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1. Introduction 

The study of neutrino oscillation in accelerator-based experiments typically reduces to the 
reconstruction of the energy of the incident neutrino which induced a reaction inside a detector, 
E v . Unfortunately, no single reaction channel completely dominates in the energy range of interest 
(roughly 0.5 < E v < 10, in GeV). Therefore, to successfully reconstruct the E v spectrum, it is usu¬ 
ally necessary to employ a combination of detailed event selection criteria and a carefully chosen 
energy estimator which reproduces the neutrino energy for a given class of interactions. But despite 
experimenters’ best efforts, the complex and overlapping nature of the reaction mechanisms usu¬ 
ally renders it impossible to avoid significant contamination from background processes. Thus, it 
is essential to have accurate predictions for both the rate and spectral shape of not only the intended 
signal process, but all probable background reactions as well. 

Historically, the available models have been constructed from admittedly crude pictures of the 
heavy nuclei from which modern detectors are constructed. Currently available simulation tools 
have performed adequately when compared to data for simple targets, such as deuterium[l]; but 
studies on heavier materials such as hydrocarbons[2] have yielded results that are murky at best. 
As a result, model uncertainty has slowly crept to the head of lists of systematic uncertainties for 
these results. Dedicated studies of the various neutrino reactions at neutrino energies of ~ 1 — 10 
GeV have thus become an indispensible part of the neutrino oscillation program. 

The MINERvA experiment, which uses the intense NuMI muon neutrino beam at Fermilab, 
is dedicated to measuring neutrino cross-sections as a function of various kinematic parameters 
on a range of nuclear targets. This information can be compared directly to models and used 
to improve or select between them. The detector itself consists of a central sampling scintillator 
region, built from strips of fluoror-doped scintillator glued into sheets, then stacked transverse to the 
beam axis; both barrel-style and downstream longitudinal electromagnetic and hadronic sampling 
calorimeters; and a collection of upstream passive targets of lead, iron, graphite, water, and liquid 
helium. Sitting as it does just upstream of the MINOS near detector in NuMI, MINERvA is also 
able to make use of the former as a muon spectrometer for most muons which exit the back of 
MINERvA. Much more detail on the detector design and performance is available elsewhere.[3] 

2. Initial-state models: CCQE 

Both of the event generators in widespread use in by modern neutrino experiments, GENIE[4] 
and NEUT[5], begin from a model which pictures the nucleons within a nucleus as non-interacting 
fermions in a relativistic gas (the Relativistic Fermi Gas model, or RFG), and select post facto 
from the plethora of available attempts at modeling interactions between them. The extent to which 
these approximations are valid can be examined directly by making measurements of the simplest 
charged interaction neutrinos experience, that of the charged-current quasi-elastic (CCQE) process, 
and comparing the data to predictions from the various models. In CCQE scattering of muon 
neutrinos, a W boson is exchanged between the neutrino and a neutron, resulting in a negatively 
charged muon and a proton in the final state: v^n —> p p. (The analogous process for muon 
antineutrino scattering reverses the lepton number and isospin: v^p —> p ! n.) We analyze this 
reaction by selecting candidate events in two different ways. 
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The first analysis[6, 7] requires an event with a reconstructed track exiting MINERvA and 
matching to a track in the MINOS near detector: this is the muon candidate. (We restrict the 
sample to events which are analyzed as having the correct charge appropriate to neutrinos of the 
helicity selected by the NuMI beam configuration; samples were collected in both neutrino- and 
antineutrino-enhanced beams.) We furthermore insist that there be no more than two isolated show¬ 
ers away from the muon vertex in the event for neutrino mode (only a single shower for antineutri¬ 
nos). With the muon energy and angle thus obtained, we are able to apply the following well-known 
estimators for the neutrino energy E v and square of the four-momentum transferred to the nucleus 
Q 2 , assuming that an event was a CCQE interaction on a stationary nucleon (with m pnp the various 
masses, and E /, a semi-empirical binding energy which accounts for the nuclear potential): 


ml - (m p - E h ) 2 - mf + 2 (m p - E b E^) 

(2.1) 

2 (m p -Eb-E^+p^ cos ) 

Qqe = 2e v E {Eg ~ Pg cos ) - ml 

(2.2) 


Our last selection requirement is an upper limit on the total non-muon (“recoil”) energy outside 
of a region immediately surrounding the muon vertex (for more on which see further below), which 
we make as a monotonically increasing function of the reconstructed Q 2 . Once reconstructed, we 
normalize the predicted backgrounds to the data by fitting them in a sideband region in the recoil 
energy. Then, we subtract the predicted background counts from the data, un-smear the muon 
energy and angle using an iterative Bayesian unfolding technique. Finally, we can normalize by 
the predicted efficiency and flux, and the number of neutrons (protons for anti-neutrinos) in our 
fiducial volume region to obtain a cross-section. Eq. 2.3 gives this in bins i of Q 2 , using e for 
efficiency, for the flux (integrated across the energy range of interest), T n for the number of 
targets, A,- for the bin width, and U/j for the predicted migration from bin i into bin j: 
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In practice, we often work with a shape-only version of the cross-section, which we can compare 
to the predictions of various models without the additional uncertainty of the flux prediction influ¬ 
encing the comparison overmuch, as in fig. 1 . The model most strongly preferred by the data in 
both v and v cases is not GENIE’s simple RFG alone, either with the world average of the free 
parameter M& (0.99) or with a larger value preferred by fits by MiniBooNE (1.35). Instead, the 
data best favor an RFG that is modified using empirical corrections to the cross-section derived 
from measurements of electron-nucleus scattering (which are thought to be necessitated by inter¬ 
nucleon correlations within the nucleus). The substantial difference from the RFG implemented in 
the GENIE generator (15-20% in some regions) underscores the necessity of moving beyond the 
unadorned RFG model. Collaborators on the T2K experiment have begun using these data to in¬ 
form the choice of models used in the NEUT generator to simulate neutrino interactions as well. [8] 


To further investigate the hypothesis of correlated nucleon initial states, we also examine the 
recoil energy spectrum in the vicinity of the reconstructed muon vertex (with displacement r < 
300mm for v and r < 1 00mm for v) in candidate signal events. Fig. 2a illustrates the discrepancy 
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Figure 1: Ratios of the differential cross-section vs. Qj L for MINERvA data and several models 
to the usual RFG model as implemented in GENIE. Other models considered: an alternative RFG 
implementation from the NuWro generator with different values of the free parameter A/ 4 , a spec¬ 
tral function (SF) model, and the RFG coupled with transverse enhancement (TEM; described in 
the text). The figure shows antineutrino scattering; neutrino scattering is similar. 


between the data distribution and that predicted by the GENIE RFG model for neutrino scattering; 
it also shows the improvement garnered by adding an additional proton’s energy to 25% of the 
interactions in the simulated sample. (We performed a fitting procedure between the data and 
simulation to determine the fraction of events that should be modified; the result was (25 ± 10)%.) 
We note the constrasting behavior of anti-neutrino scattering, shown in fig. 2b: here the discrepancy 
is marginal, and the fitting procedure favored a modification of (—10±7)% of events, which is 
consistent with no modification at all. Taken together, these results lend credence to the picture of 
initial-state correlations between neutron-proton pairs in some events: neutrino scattering converts 
a neutron to a proton, resulting in a post-reaction pp pair, and thus an excess proton which we 
observe in our data but not our simulation; antineutrino scattering, on the other hand, creates a nn 
pair, and the additional neutron typically (but not always) goes undetected in MINERvA, so the 
distributions agree better without modification. 

Our second CCQE analysis seeks to approach the kinematic measurement from a different 
direction. Most of the selection criteria are analogous to those detailed above. But instead of 
demanding that the muon track be analyzed in MINOS, this analysis simply requires the muon 
candidate to exit the detector (itself a reasonably pure selector for muons). More significantly, for 
this measurement we also require that the final-state nucleon’s trajectory be reconstructed and that 
its dE jdx profile positively identify it. (Since neutrons are typically unobserved in MINERvA, 
this therefore restricts the measurement to neutrino mode only, where the final-state nucleon is a 
proton.) We then can make the comparison corresponding to fig. 1, except that we choose to use 
the proton ’’s kinematics (including its kinetic energy, T p ) to reconstruct Qq E : 

Qqe, p = (M„ ~ Eb? + 2 (M„ - E b ) (T p - M n - E b ) - M 2 (2.4) 
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Figure 2: Comparison of GENIE RFG prediction (before and after fitting) of energy near recon¬ 
structed CCQE event vertex vs. MINERvA data. The fit used for “Best Fit” adds an additional 
proton’s energy and is further described in the text. 


Fig. 3 shows this comparison. Here, intriguingly, most of the models fit the data similarly well, 
with the GENIE RFG model being the slight favorite. No single explanation for the stark contrast 
with the result in the muon kinematics (fig. 1) has yet gained favor, though suggestions include the 
potential inaccuracy of GENIE’s background model as well as a need for better-developed models 
of multinucleon effects implemented in generators. [9] 
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Figure 3: Ratios of the differential cross-section vs. Q 2 for MINERvA data and several models 
to the usual RFG model as implemented in GENIE. “LFG” means “local Fermi gas”; “RPA” is a 
random phase approximation!! 1]; “Nieves” is the nucleon correlation model of Nieves et al.[12]; 
TEM is a transverse enhancement model [13]. 


Work on additional cross-section results for the CCQE process is currently underway. A dou- 
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bly differential cross-section (jjppjy ) analysis, which will examine the validity of the model pre¬ 
dictions in more detail, is in preparation. A preliminary result for the electron neutrino CCQE 
cross-section (total cross-section vs. energy, as well as differential cross-section vs. 0 e , E e , and Q 2 ) 
was also presented at this conference.[10] 

3. Final-state models: 71^ production 

Since charged pions are both light and strongly interacting, single charged pion production by 
neutrinos (in MINERvA, V^CH —? /i n + X, where X is anything but a n + ) is an ideal reaction 
channel for the study of the final-state (re-)interactions (FSI) of hadrons after a neutrino interac¬ 
tion. Comparisons of detailed measurements from the MiniBooNE experiment to the predictions 
of various generator models have raised questions about whether current approaches to FSI mod¬ 
eling are even valid.[14] Thus, further measurements are extremely desirable for confirmation or 
disambiguation of this situation. 

For our analysis[15], we select a sample of events containing a charged-current interaction 
with a single charged pion in the final state by requiring a MINOS-matched muon candidate and 
one or two other reconstructed tracks. Exactly one of these tracks must both be consistent with a 
pion’s dE/dx energy loss profile and also be associated with spatially close, time-delayed activity 
consistent with a decay electron. We also cut on the reconstructed hadronic invariant mass W = 
\JMj, — Q 2 + 2M p Ehad (assuming target nucleon at rest) of the reaction, requiring W < 1.4 GeV, so 
as to select a sample purer in single pion events (typically resulting from decay of the A resonances 
with W ~ 1.2 GeV). 

After extracting shape-only cross-sections vs. pion kinetic energy after the manner of eq. 2.3, 
we are able to compare our data to models which use various different FSI approaches or tunings. 
This is shown in fig. 4. Models with a full treatment of FSI are strongly preferred by our data, 
as opposed to models with no or partial FSI simulation. Once again, collaborators on the T2K 
experiment are using these data to constrain the FSI parameters used in the NEUT generator to 
simulate neutrino interactions. 

4. Cross-section models: coherent n- 

Finally, we study coherent pion production, in which the momentum transferred from the 
neutrino is entirely used to create a new particle, and the struck nucleus remains in its ground 
state. In its charged form (v u 4 —> ji n + A, or, exchanging leptons for antileptons in antineutrino 
scattering), this process is a background for v^ CCQE processes; the neutral form (v u 4 —y v Lt 7i°A) 
can mimic v e CCQE when one of the decay photons from the jc° is unobserved. 

For the MINERvA analysis[16], we begin selecting events consistent with charged-current 
coherent scattering by demanding they have exactly two reconstructed hacks. We further refine the 
sample by cutting on the reconstructed momentum transferred to the nucleus, |f|, and the energy 
observed within a “box” of side length 40cm and depth 15cm centered around the reconstructed 
event vertex. As before, we construct cross-sections in pion angle according to eq. 2.3, and then 
compare them to the most commonly used model for coherent pion production, that of of Rein 
and Sehgal, as implemented in GENIE[4]. This is given in fig. 5. We observe that the measured 
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Figure 4: Comparison of various predictions for do/dT n with and without FSI vs. MINERvA 
data. GENIE FSI. NEUT, NuWro, and GiBUU have full FSI models; GENIE No FSI and ACS 
have no or partial FSI. (Models are described in more detail in ref. [15].) 

distributions are significantly more forward than the model in both neutrino (shown) and anti¬ 
neutrino scattering (not shown here). Further discussion of this result and its implications were 
given in a separate presentation at this conference. [17] 



0^ w/r to Beam (Degrees) 


Figure 5: Comparison of the Rein and Sehgal model implemented in GENIE for do/dO K vs. 
MINERvA data. 


5. Conclusions 

Measurements from MINERvA provide valuable information in the quest to improve nrod- 
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els of neutrino interaction processes and neutrino energy estimation. We find inadequacies in the 
simple RFG picture usually used for CCQE, as well as in the underlying cross-section model for 
coherent pion production, while models of FSI currently used in generators appear to be sufficient 
to model our data. These results are already being applied by the T2K collaboration in their gen¬ 
erator model tuning, and we anticipate that other generator efforts can benefit from similar studies. 
We also expect both forthcoming and planned work (in the low-energy beam configuration used 
in results discussed here, as well as the medium-energy tune in which MINERvA is currently col¬ 
lecting data) to continue to improve the fidelity of neutrino interaction modeling, and therefore the 
quality of reconstructed energy distributions. 
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